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Abstract 



We present results for the Bs meson decay constant /b^ from simulations 
at three lattice spacings in the range a^^= l.lGeV to 2.6GeV using NRQCD 
heavy quarks and clover light quarks in the quenched approximation. We 
study scaling of this quantity and check consistency between mesons decaying 
from rest and from a state with nonzero spatial momentum. Cancellation of 
power law contributions that arise in the NRQCD formulation of heavy-light 
currents is discussed. On the coarsest lattice the Dg meson decay constant 
fos is calculated. Our best values for the decay constants are given by Jb^ = 
187(4)(4)(11)(2)(7)(6) MeV and fo, = 223(6)(31)(38)(23)(9)(t?) MeV. 
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I. INTRODUCTION 



The B meson decay constant, fs, defined through the hadronic matrix element 

{0\A,\B,p)=p,fB, (1) 

is one of the phenomenologically relevant quantities whose determination has relied heavily 
on lattice calculations. Together with the bag-parameter Bb, the decay constant is an 
important ingredient in analyses of B^B^ mixing phenomena and studies of CP violation in 
the neutral B system. Several quenched lattice determinations of Jb have been completed in 
recent years using different approaches to simulating the heavy 6-quark on the lattice 
The first dynamical results are also starting to appear P-pT|,|^]. The most extensive quenched 
studies have been carried out using either the non-relativistic (NRQCD) formulation |12| or 



the heavy clover approach |[T^ to 6-quarks. Results for fs and J'b^ from these two methods 
are in good agreement with each other. 

In this article we investigate once again the NRQCD approach to leptonic decays of heavy 
mesons. We wish to gain further insight into this method, which will strengthen confidence 
in our previous results and also allow us to assess its potential when going on to other 
important areas of heavy quark physics such as semi-leptonic decays. Our analysis is based 
on quenched results for the Bs meson decay constant from simulations at three lattice 
spacings spanning = l.lGeV - 2.6GeV. We use NRQCD 6-quarks, clover light quarks 
and the Wilson plaquette gauge action. Simulation parameters are summarized in Table |. 
The (3 = 6.0 Jb^ results in sections III & IV are taken from reference 0. The configurations 
at /? = 5.7 and /? = 6.2 plus the light propagators at /3 = 6.2 were provided by the UKQCD 
collaboration. The (3 = 6.0 data presented in section VI also uses UKQCD configurations and 
light propagators. Some preliminary results at earlier stages of this project have appeared 
in 



mpi. 



Our first goal was to check for independence of our final fs^ value from the lattice spacing 
used. This is particularly important here, since we employ an effective theory, NRQCD, for 
heavy quarks which precludes extrapolations to zero lattice spacing, a —>■ 0. There are many 
advantages to NRQCD, including being able to simulate quarks with masses much larger 
than the cutoff {aMq S> 1) without introducing large discretization errors and the efficiency 
with which heavy quark propagators can be computed numerically. However, the theory is 
not renormalizable and furthermore aMg cannot become too small. In order for the approach 
to work, one must have a window in which the lattice spacing is large enough so that aMg 
is of 0(1) or greater, but also small enough so that all discretization errors, including those 
coming from the glue and light quark sectors, are under control. Physical results must 
already be independent of the lattice spacing within acceptable and quantifiable systematic 
errors. The scaling studies presented in this article demonstrate that such a window does 
indeed exist between = 1.1 ~ 2.6 GeV for simulations with NRQCD-clover quarks and 
Wilson glue. The systematic errors are ~ 10% in the upper half of the window and increase 
to ~ 20% towards the coarse end. By further improving the quark and gauge actions one 
should be able to improve the errors at the coarse end and even extend the window down 
below IGeV. Recently the JLQCD/Hiroshima collaboration M studied scaling of fs and 
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with a lattice action very similar to the one used in this article. Our results are in good 
agreement with theirs. 

Another feature of the use of an effective theory, is that the matrix element of current 
operators in the full theory (continuum QCD) such as the LHS of eq.(|l|) must be built up 
out of several current matrix elements evaluated in the effective theory, 

{A,)QCD = j:C^i4'^)eff- (2) 
j 

With NRQCD heavy quarks the currents are of the form, 

where the 0|/^ are operators of dimension dj = 3 + Ij {Ij > 0). Quantum corrections to the 
currents and the mixing of higher dimension operators back onto lower dimension operators 
will induce "power law" terms, i.e. terms that go as ag/ {aMq)^^^ (plus higher orders). 
Since the LHS of eq.(0) should have no power law contributions, any consistent framework 
for evaluating (</^-'^)e// and Cj must arrange for power law terms to cancel to the order 
one is working in. The second goal of this article was to investigate and quantify this 
cancellation at 0{as/ {aMq)). In our simulations we include all relevant current operators 
through dimension 4. The matching coefficients Cj have been calculated perturbatively 
through 0{as/MQ) and 0{aas) pOf. Matrix elements of the dimension 4 current operators 
will contain in addition to Aqcd/Mq physical contributions also power law terms of which 
the dominant 0{as/{a Mq)) terms are cancelled by our matching coefficients. Remaining 
uncancelled power law terms start at 0{al/ {a Mq)) and their effects are taken into account 
in the systematic errors that we quote. Our analysis, based on the temporal component of the 
axial vector current, Aq, demonstrates that, although power law contributions are substantial 
there are no problems with cancelling them and there are no delicate fine tunings involved. 
Furthermore, we find that after cancellation of 0{as/{aMQ)) terms the contributions from 
1/Mq current corrections are significantly reduced. 

The third goal of the present study was to carry out other consistency tests in our lattice 
evaluation of /^^ . The definition eq.(|l|) of the decay constant tells us that one should be 
able to obtain the same Jb^ from Bg mesons at rest and from those with nonzero spatial 
momentum up to lattice artifacts. We have simulated leptonic decays of Bg mesons with 
nonzero momenta to check this. Another project, which we will report on in a separate pub- 
lication, is determining Jb^ from the spatial component Ak, providing a further covariance 
test of calculations with NRQCD heavy quarks 



In the next section we introduce the lattice actions and current operators used in our sim- 
ulations. Section III presents results for /^^ at the three lattice spacings and scaling plots. 
In section IV we consider cancellation of power law contributions and also discuss the slope 
of fpsy/Mps versus 1/Mq as one leaves the static limit. On our coarsest lattice we are able 
to simulate NRQCD charm quarks and in section V we present results for the Dg meson 
decay constant f^^ at this one value for the lattice spacing. We describe /^^ extracted at 
nonzero momenta in section VI and section VII gives a brief summary of this article. In an 
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Appendix we provide a Table of heavy-light current matching coefficients covering a wider 
range of heavy quark masses than in [^] . 



II. THE LATTICE ACTION AND CURRENT OPERATORS 



We used the NRQCD heavy quark action density given by, 




where V denotes a two-component Pauli spinor, Hq is the nonrelativistic kinetic energy 
operator, 

A(2) 

and 5H includes relativistic and finite-lattice-spacing corrections, 

(API)' g'Al" a(AP')' 

8(A/„)3 ^ 24A/„ '''l6n(M„)2' ^ ' 

In addition to all terms we have also included the leading relativistic correction 

as well as the discretization corrections appearing at the same order in the momentum ex- 
pansion. For all three P values the NRQCD action was tadpole- improved using the plaquette 



definition of mq [H] and the Cj were set equal to unity. Table II lists the bare heavy quark 
masses aMo at which simulations were carried out together with the corresponding values 
for the stability parameter n. 

For the light quarks we used the tree-level tadpole improved clover action {csw = Wq^) 
at /5 = 5.7 and f3 = 6.0. However at /3 = 6.2 the nonperturbatively determined value 
csw = 1-61 was employed The k values used in the simulations are listed in Table I. 



At /3 = 5.7 and /? = 6.2 we employed a single k close to Ks fixed through the K meson. At 
f3 = 6.0 data was obtained at three light k's, 0.1369, 0.1375 and 0.13808, and results were 
then interpolated to Ks 0] . In Table I. we also list the lattice spacing determined from the p 
mass and the bare dimensionless heavy quark mass aM^ corresponding to the b quark. The 
latter was determined from a combination of perturbation theory for mass renormalization 
and simulation results for binding energies. These procedures for fixing aM^ have been 
explained in several previous publications |3,15]. The simulation parameters given in this 



paragraph apply to results presented in sections III, IV and V. The (3 = 6.0 simulations 
described in section VI use different parameters (see section VI and Table VII). 

For the temporal component of the heavy-light axial vector current eq.(E|) becomes, 
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with 



{^o)QCD = T.C\^'\ji'^) + 0(1/MJ), (7) 

j=0 



= (ll5loQ, 



(8) 



The heavy quark 4-spinor Q has the NRQCD 2-spinor as the upper two components, 

zero for the lower two components and obeys 70Q = Q- Jq'^ and Jq contribute already at 

(2) 

tree-level whereas Jq appears only at one-loop. 

Higher dimension operators appear in lattice effective theories when one, for instance, needs 
to consider relativistic corrections. They also occur whenever lattice operators are improved. 
Improvements to current operators take on the form, 

jU)disc ^ a^jO^\ (9) 

where again, as in eq.(|^), the O^^^ are operators of dimension dj = 3-|-/j. In the present case 
of NRQCD-clover currents, one has the identical set of possible dimension 4 operators O^^^ 
in eq.(^ and eq.(|^). Hence a one-loop mixing and matching calculation to full continuum 
QCD involving the above J^^ will not only tell us how the I/Mq currents contribute but 
will also identify the discretization corrections that come in at one-loop. In other words a 
consistent matching through 0{as/ Mq) should automatically be consistent through 0{aas)- 
In reference it was found that for each /i only one of the J^-?)'^*^^ was relevant at one-loop. 
Matching to continuum QCD could be achieved by improving J^°^ in the following way, 

JT Ji^^+C^sCA,jjf-\ (10) 

with 

4"'"^ = -aqj-^lol5l,Q = 2aMoJj^\ (11) 

and (ao = 1-03. The effect of Jl^'""^ can be taken into account by expressing it in terms of 
jj^^ and absorbing as2aMo(Afj, into the coefficient 6*2'^^''. 

At P = 5.7 and /? = 6.0 1/Mq current corrections that appear at tree- level have also been 
studied. For Aq these are, 

J^'^ = ^g757oD^g, 



-/o'^ = 4^^757oa-EQ, (12) 
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where ex. = 707 and S = diag{(T, cr). We will comment on matrix elements of these currents 
and their Mq dependence later in section V. However, since the 0{as/ Mq) full matching cal- 
culation has not been carried out yet, 0{as/ {a MqY) and 0{as{I^QCD I Mq) / {a Mq)) power 
law contaminations contained in these matrix elements will not be cancelled and we prefer 
not to include them in our scaling analysis of fs^ ■ In reference tree-level contributions 
from (|l^) were included. The total effect was a ~ 3% decrease in fs^ , well within the 
systematic errors quoted there for uncancelled power law terms. In the (3 = 6.0 results given 
in the present article, however, we have removed the 1/Mq tree-level current effects. 

Based on the lattice actions employed, the improvement of J^^ and the consistent matching 
to continuum full QCD through 0{as/ Mq) and 0{aas) in the currents, one can now list 
the expected remaining dominant systematic errors in our evaluation of (A^)qcd- These 
are collected in Table III and their size, in percent, estimated for each (5. In making these 
estimates we use Kqcd ~ 400Me1/, = 5.7) ~ 0.33, = 6.0) ~ 0.24, = 6.2) ~ 
0.20 and approximate aMq with oMq of Table I. The as values are close to ap{q* = 1/a), 
where ap is the coupling introduced in reference closely related to «y which is based on 
the heavy quark potential [^,^. Discretization errors come in at 0((aAQC7£))^) from the 
currents and at /5 = 5.7 and P = 6.0 also at 0{as qAqcd) from the light quark action. These 
latter errors arise from using a tadpole-improved csw rather than a non-perturbative csw 
(or a Csw value fully corrected at one-loop) and are therefore proportional to the difference 
of these values for csw- Such differences and, hence, also the O^asaAqcD) errors are small 
at /? = 6.0. Relativistic errors are at O^^Aqcd/Mq)"^) from the currents. There are also 
0{as Aqcd/Mq) errors from using the tadpole-improved coefficient for the a ■ B term in the 
NRQCD action instead of the full 0{as) coefficient matched to continuum QCD. This error 
is proportional to the l/Mg slope induced in fs by the a ■ B term. The investigations in 
found this slope to be small, approximately 0.25* Aq^d/^q- We fold in perturbative errors 
separately for O(a^) and 0{a1/ {aMq)) since these have different /? dependence and also to 
indicate that perturbation theory will eventually break down for small aMq, i.e. for large 
enough (3. We note that in the absence of tree-level 0{l/{Mq) current matrix elements, 
there are no 0{as/ {aMqY) errors in our calculations. As long as the same action is used 
in both the simulations and the perturbative matchings, all 0{as/ {aMqY) terms cancel 
between matrix elements and matching coefficients. Table III also lists other systematic 
errors (uncertainty in Kg and fixing of from the p mass) and the statistical errors in the 
determinations. Not listed are uncertainties due to quenching. Recent studies indicate 
that these last corrections can be substantial, at the 10 ~ 25% level ||9HlT . 



III. SIMULATION RESULTS AND SCALING OF fs. 

In this section we describe /^^ calculated from Aq and investigate scaling of this quan- 
tity. We start with simulation results for current matrix elements in lattice units before 
renormalization. Table |V| summarizes our results at /? = 5.7 and (3 = 6.2 for 

3/2 

a^/^ Z^^) = (0| \PS). (13) 

Mps 
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Data are presented for several heavy quark masses with the hght quark mass fixed to Kg of 
Table I. The corresponding simulation data for f3 = 6.0 can be found in The matrix 
elements were evaluated in the meson rest frame. ( Jq^^ ) is identical to ( Jq^^ ) in this frame 
and hence is not listed separately in the Table. For /3 = 5.7 we also list results for the higher 
order currents of eq. ([12|) . 

Results in physical units {GeV^^'^) for 



fps = E (14) 

j=0 

are summarized in Table V, at tree- level {Cq^"^ = Ci^^'^ = 1, C2^°^ = 0) and with one- loop 
matching for three different scales to fix = ap{q*). At /3 = 5.7 and (3 = 6.0 we can directly 
take over the matching coefficients from Appendix A, which assume tadpole improvement 
in both the NRQCD and clover actions. At /5 = 6.2 where the non-perturbative csw is used, 
some slight modifications are necessary. We first remove the tadpole improvement term 
coming from the light quark wave function renormalization (denoted C'^^ in [Q) from the 
matching coefficients. We then use k rather than k = k x uq to rescale the light quark field 
in the currents. As a result the entries in Table IV and the tree-level results in Table V are 
enhanced for f3 = 6.2 by approximately a factor of Uq This enhancement is compensated 
for to a large extent upon going to renormalized matrix elements (the latter three columns 
in Table V) due to larger perturbative matching coefficients. Alternatively, we could have 
used another light quark action at j3 = 6.2 which has the first derivative terms divided 
by Uo and a clover term with coefficient C^^^^^^' /uq. This would have allowed us to use 
tadpole- improved perturbation theory even for the light quark propagator, similar to what 
was done at the other (3 values. We tested this option and found that the final answer for 
the decay constant changed by 4 ~ 5%. This difference is, as expected, of the same order 
as other O(a^) corrections not included in the present analysis. 

The scale q* has not been calculated yet for the matching coefficients with NRQCD heavy 
and clover light quarks, q* is known for static heavy quarks combined with massless Wilson 
light quarks ^6| and also for light-light currents |2^ again using massless Wilson fermions. 



These authors find q* = 2.18/a and q* = 2.3 /a respectively, suggesting very mild dependence 



of q* on aMq. Results also exist for clover light-light currents p8[ showing some reduction in 
q* relative to the Wilson results for these currents. In the absence of an explicit calculation 
for the action and currents employed in the present simulations, we use the q* = 2/a numbers 
to quote central values and arrange for perturbative errors to cover the spread in results due 
to uncertainty in q*. Table VI shows our final numbers for fs^ interpolated to the physical 
heavy quark mass Mq for the three lattice spacings. The Particle Data Group's 5° mass, 
5.369GeV, was used to convert from fps\/Mps to /b^. The errors follow those listed in Table 
III and correspond from left to right to statistical (plus interpolation in Mq), discretization, 
perturbative, relativistic. Kg determination and errors in a~^{m 



pi 



After adding all errors in quadrature, we plot /^^ versus a in Fig. 1 and compare with 
results from H]. One finds good scaling within errors and excellent agreement between the 
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two collaborations. These results also agree well with quenched fs^ determinations using 
other heavy fermion formulations on the lattice [|l|. 

Table VI and Fig.l represent one of the main results of this article. Working around of 
~ 2GeV with NRQCD-clover quarks and Wilson glue and matching through 0{as/ Mq) and 
0{aas)i one can obtain reliable quenched numbers for B meson decay constants with ~ 10% 
accuracy. To do better at these lattice spacings one would need to go to higher orders in the 
matching procedure. At our coarsest lattice spacing the systematic uncertainties are larger. 
Going to more highly improved actions would help here. Nevertheless, the /? = 5.7 results 
are consistent with our best numbers. Simulations with ~ IGeV can already give a 
good indication of continuum physics. 

Since we cannot extrapolate to a — 0, in order to quote our best estimate for quenched 
there are several choices. One approach would be to average over the [3 = 6.0 and 
f3 = 6.2 results. Alternatively one could, after having verified scaling behaviour, pick the 
point with the smallest statistical and systematic errors. Based on Table III this corresponds 
to the (3 = 6.2 result. We adopt the second method here and present as our final best value 
for the Bs meson decay constant in the quenched approximation, 

/b, = 187(4)(4)(ll)(2)(7)(6)MeV, (15) 

with errors as described in Table VI. 

If the value of /^^ at the physical 6— quark mass were our only goal then we would be done at 
this point. However, if further details are of interest, such as " the true size of contributions 
from 1/Mq current corrections to fs " or " the slope of fpsy/Mps versus 1/Mq ", then one 
needs to scrutinize more carefully the contributions from individual current matrix elements 
( Jq'^ ) and matching coefficients Cj to the RHS of eq.(|^). We turn to these issues in the 
next section. 



IV. POWER LAW TERMS AND Kqcd/Mq CORRECTIONS 

In this section we investigate further the properties of matrix elements of the 1/Mq 

current corrections. We focus on the current Jq^"* = 2l|;^757o 7"^ This operator is 
introduced in order to take into account p/Mq effects which we know are present in the full 
theory, i.e. in the LHS of eq.(^. Here "p" is the momentum of the quarks inside the meson 
and should be of 0{Aqcd)- Once introduced into the effective theory, however, this current 
will mix with other currents in the theory. Its short distance effect will be to renormalize 
the zeroth order current and one expects the matrix element ( Jg^^ ) to develop a component 
proportional to {Jq^^)- On dimensional grounds the proportionality coefficient will be of 
the form const/ {a Mq) and one has, 

( Jq^^ ) = — —{ Jn^^ ) H 9£R term + higher orders. (16) 

uMq Mq 

The first term on the RHS of ([T6| ) is an example of a power law contribution to a matrix 
element. They are unavoidable in a quantum effective field theory where one has mixing 
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between operators of different dimensions. However, the matching procedures between the 
full and the effective theories should remove these artifacts. Since we do the matching 
perturbatively, power law contributions can only be subtracted to a given order in as- 
As mentioned several times already, in our present calculations we cancel terms through 
0{as/ {a Mq)). To see precisely how this cancellation comes about we need to dissect the 
matching coefficient Cj=Q. Using the notation of ||20| one can rewrite eq.(^ as. 



with 



^0)qCD = (1 +«sP0) (4°^) + (1 + «.Pl) (J^'^) +«sP2 (Jf ), (17) 



po = 5o - ^(C, + Cq) - Coo - Clo = Po - Clo. (18) 



Po is the one-loop contribution to Cq and consists of parts coming from continuum perturba- 
tion theory (-Bo)? from light and heavy quark wave function renormalizations on the lattice 
[Cq and Cq) and from mixing between the currents Jq"'^ (Coo and Cio)- Coo is the feedback of 
Jq°'* onto itself and Cio incapsulates the projection of Jg^'* back onto Jq^\ precisely what one 
is looking for in the first term of ([l6|). Picking out the relevant terms, one can write 



(4°') + ^*-™+ ... (19) 



;i + CXs Po) (4°^) + (4'^) = (1 + «s Po) ( jf) 

Clo /. 



+ 
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aMQ 

(0) 



The power law contribution — as Cio ( -^o"'' ) has now been reduced from an 

0{as/ {a Mq)) to an 0{a1/ {a Mq)) effect through the a^Cio matching term. It is conve- 



nient to define a subtracted Jq 



(1) 

(4'^)s„5=(J,S'^)-«sCl0(^r), (20) 

which is equal to the last two terms in ([T9|) up to higher order corrections. It is ( Jq^"* )sub 
rather than ( j'^^ ) that should be used to assess the importance of 1/Mq current corrections 
to fs,. In Fig.2 we plot the absolute values of f^^^VM, f^^\fM and /^VM (related to 
(Jq^'*), {J'q^) and {JQ'^)sub in the usual way) versus 1/Mps at = 6.0. One sees that 
|/j^|,|A/M is considerably smaller than \f^^^\/M and amounts to only ~ 4% to ~ 5% of the 
lowest order f^^^y/M term. This is of the order of systematic errors such as 0(af) errors 
remaining in our calculation, so we cannot give a more precise estimate of the true Aqcd/Mq 
contributions coming from 1/Mq currents. Similar results are obtained at the other lattice 
spacings. A corresponding plot similar to Fig.2 for f3 = 5.7 will be discussed in the next 
section. At all three (3 values the power law contributions make up 50% to 80% of ( Jg^^ ), 
depending on q*, and once they are subtracted away one is left with only a small remainder. 
We note that the subtraction in (pO|) took place automatically in the Jb^ calculations of the 
previous section. The fact that /j^),v^M turned out to be small did not pose any problems 
for the power law cancellation itself at O^as/iaMg)). 

Another quantity where having information on Aqcd/Mq corrections is important, is the 
slope of fpsVMps versus 1/Mq. Writing, 
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fps^M^s = (fps^/M^) [1 + ^ + 0(1/M2)], (21) 

\ / Stat Mq 

one is interested in the slope A^^. The bulk of the effect will come from 1/Mq terms in the 
NRQCD action which are responsible for the Mq dependence of the zeroth order current 
matrix element ( Jq'^ ), and from the one- loop correction po to this matrix element. We have 
just seen that contributions from l/Mg current corrections are small. Nevertheless, when 
simulation data are fit to the form (|2T|) to extract A^^ one might worry about contaminations 
due to uncancelled power law terms of 0{a1/ [a Mq)). If A^; is the result of a fit then one 
can estimate a percentage error of a1/{aKsi) in this slope. Putting in some numbers this 
leads to A^^ceF] ^ uncertainties depending on (3, where A^^ must be inserted in GeV's. If the 
slope is less than Aqcd (0.4 GeV) these errors can become quite large. 

There is, however, another reason why it is difficult to extract an accurate A^; from our 
simulations. In order to calculate the slope we need to look at data at masses larger than 
the b quark mass. The one-loop terms in the matching coefficients increase as one goes 
towards the static limit and as a consequence results for fpsV^PS become much more 
sensitive to the value of as i.e. to q*. This can be seen in Fig. 3 where fpsV^PS is plotted 
versus 1/Mps for q* = 1/a and q* = Tc/a. One sees that there is considerable dependence of 
the slope on q*. We will quote a central value for the slope based on the /3 = 6.0 data using 
q* = 2/a and allow for a change in q* in our error. We note that in order for the expansion 



in ([21| ) to be valid there should be only power dependence and no logarithmic dependence 
on Mq. We satisfy this condition by fixing the logarithms in the matching coefficients to 
log(aMo^). Our best estimate for A,; is then: A,, = -0.67(9)(;3l)(12)GeV. The first error is 
the statistical fit error to the q* = 2/a data and the second comes from the upper limit of 
errors on the magnitude of the slope from fits to q* = ir/a data and the lower limit of errors 
in fits to g* = 1/a data. The third error is an estimate for and a'^/{aMQ) corrections. In 
Fig. 4 we plot the P = 6.0 results together with data from other /3 values all using q* = 2/a. 
One sees that slopes extracted at the two finer lattice spacings are consistent within our large 
errors. Results from j3 = 5.7 lie mostly parallel to the (3 = 6.0 and (3 = 6.2 data, however 
statistical errors are too large near the static limit to enable a useful estimate for the slope. 
Around the physical B meson the f3 = 5.7 data points lie ~15% high, in agreement with 
results in Table VI and consistent with our estimates of systematic errors in Table III. 

In contrast to J'b^ which we demonstrated can be calculated with ~ 10% errors, we find that 
specific 1/Mq corrections are harder to determine with comparable percentage accuracy. For 
the two quantities discussed in this section, the 1/Mq corrections coming from the current 
corrections and the slope of fpsy/Mps as one leaves the static limit, the main reason is 
uncertainty coming from higher order perturbative corrections. 



V. SIMULATIONS AT THE CHARM QUARK AND fo. 

In this section we discuss the behaviour of the decay constant and the current matrix 
elements for heavy quark masses in the charm region. For the coarsest of our lattice spacings 
the charm quark can be reached at a bare mass of aMo^ = 0.87(6) (3) (ijg) pf. 
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The size of the lattice current matrix elements for different values of the heavy quark mass 
is compared in Fig. 5 . This figure is similar to Fig.2, however at /? = 5.7 we cover a much 
wider range of mass values. We also plot l/Mg current matrix elements. The figure shows 

the value of |/i^f,| to be subleading even in the charm region, however it has grown from 
4% of /(°)a/M at the Bs to 14% at the D^. For aMg ^ 1 the value of the 1/M^ current 
matrix element d/*-^-' + /''^^ + f^^^\)\/M is equal to \f^^^\\/M and the two lines in Fig. 5 cross. 
After the discussion in the previous section, this is not unexpected and provides evidence 
that this current is also dominated by power law contributions. The part of the current 
{Jq'^ + Jq'^ + Jq^^) proportional to {Jq ) diverges as as/ {aMqY . For oMq ~ 1 this results in 
a suppression of 0{as) with respect to (Jq'^), which is the same suppression factor as the one 
for ( Jq^'*) . In order to avoid uncertainties of 0{as/ (oMq)^) in the final fo^ , we do not include 
the matrix elements of 1/Mq currents. By doing so, we have to accept an uncertainty of the 
size of Aqq<j^/Mq, which is 10%, for the final result. This is a conservative estimate, however. 
Based on our experience with Aqcd/Mq current corrections, we can expect that most of 
Q j(3) _|_ y(4) _|_ plotted in Fig. 5 will eventually be cancelled by matching coefficients 

at Oias/iaMQ^). So, even at the Dg, the true Aq^j^/Mq current corrections should not be 
a large effect. The other Aq^j-,/Mq contributions in fo^ come from the Ag^^/Mg terms in 
the NRQCD action, and these are already included in the present simulations. To estimate 
the total NRQCD error we add to the 10% another 2% coming from 0{as Aqcd/Mq) errors 
in the NRQCD action (see discussion at the end of section II). 

For the perturbative uncertainties we assign 11% for a^. This is larger than the uncer- 
tainty of arising from the variation of q* in a range from 1 to tt. For the residual power 
law contributions of (9(a^/aMg) we assign an uncertainty of 13%, which leads to 17% for the 
total perturbative error, when added in quadrature. The 13% uncertainty is comparable to 
the size of the entire {JQ^^)sub matrix element. For the uncertainty in the strange quark mass 
we assign 4% as in the case of fs^, which turns out to be negligible compared to the other 
sources of uncertainty. It is interesting to note that the uncertainty arising from the scale 
a{mp) almost cancels between the effect this has on the conversion of the matrix element 
into physical units and the shift of the bare charm quark mass. 

For our final result we obtain 

/d^ = 223(6)(31)(38)(23)(9)(1?) MeV ^ /z?, = 223(55) MeV . (22) 

The uncertainties are listed in the same order as in Table |IT| with the two discretization 
errors added in quadrature. This result is in very good agreement with the lattice average 
of 220^20 MeV quoted in ||29| for quenched fo^ ? however the uncertainty we assign is larger. 



It is also interesting to compare with direct experimental determinations of fo^ ■ In a recent 
review article the authors of reference |3^ present a world averaged experimental value of 



ff)^ = (241 ± 32)MeV, which is very consistent with the lattice numbers. Recall, however, 
that the lattice results given here are all in the quenched approximation. 
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VI. fB, FROM FINITE MOMENTUM MESONS 

This section investigates the defining matrix element eq.(l) for nonzero spatial momen- 
tum p at two lattice spacings, P = 6.0 and /? = 6.2. Our goal is not to calculate fs^ more 
accurately than we have done in section III from B mesons at rest, but to test our ability 
to simulate hadrons at finite spatial momenta and obtain consistent results. At each f] we 
work at one value for the heavy quark mass and we do not attempt to extrapolate to the 
physical b quark. Furthermore ai (3 = 6.0 we use a different set of configurations from 
the previous sections. The new configurations and light propagators were provided by the 
UKQCD collaboration. In Table VII we summarise relevant action parameters. 

We accumulated data for four nonzero momenta, (0,0,1), (0,1,1), (1,1,1) and (0,0,2) in units 
of 27r/aL averaging over all equivalent momenta ±p^., ±py etc. From the difference of falloff 
energies of meson correlations with and without momenta, 

SE{p) = E,,^{p)-Es^m{0), (23) 

one can use a relativistic dispersion relation to define a mass for the heavy meson, usually 
called Mkin, 

M,,n = {p^-SE^)/{26E). (24) 

Our first test is to check the extent to which M^in is independent of p. This nonperturbatively 
determined mass can also be compared with another mass Mpert based on the perturbative 
pole mass for the b quark combined with an energy shift Eo and simulation results for 

Esimip = 0), 

Mpert = ZmMo -Eo + Es^m (0) . (25) 

We use one- loop perturbation theory for and Eq. The one- loop coefficients and q* for the 



combination (ZmMo — Eq) are given in |T5|. In Table VIII we list results for M^in and Mpert 
in lattice units. The errors in Mkin are statistical and for Mpert we give 0{al) perturbative 
errors. The statistical errors in Mpg^ coming from Egi^iO) are negligible compared to the 
perturbative errors. The quality of our signals for finite momentum mesons depended rather 
strongly on the smearings employed (see |T^ for details of smearings). In some instances. 



e.g. in the P = 6.0 data, clearer signals were obtained for the highest momentum than for 
some of the lower ones. We believe this is because the smearings used were better suited for 
that momentum. For momentum (1,1,1), on the other hand, we could not extract any useful 
results at /? = 6.0. Fig. 6 plots aM^^ versus the momentum and compares with aMp^rt- 
Mkin is reasonably independent of p for the range considered and agrees with Mpert within 
~lo". The statistical errors at some of our data points, however, are still considerable for 
Mkin and they could be hiding systematic discretization corrections. 

The next quantity of interest is the ratio of the decay matrix element for a heavy meson with 
momentum p to that of a heavy meson decaying at rest |^. More specifically we consider. 



(0|^|P5 , p)/^E{p) ^E{p) 

(26) 



{0\Ao\PS , p = 0)/y/M^ y/M^'' 
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where E{p) is the total energy of the meson. The RHS of eq.(p6D is a very slowly rising 
function of p, never getting very much above 1 for the momenta involved. In our simulations 
we take the following expression for the LHS, working to the same order in 1/Mq as in the 
previous sections. 



Whereas Mkin, discussed above, checks for the correct relativistic dispersion relation in the 
energies, R{p) tests the p dependence of amplitudes. In Table IX we give results for R{p). 
We also list R'^^\p), the analogous ratio using only the zeroth order current Jq^\ 



One sees only a small change at less than the la level between R{p) and R^^\p). Most of the 
effect of one-loop matching and 1/Mq current corrections is canceled in this ratio. In Fig. 7 
we plot R{p) versus the momentum. The full line is the expected behaviour coming from 
the RHS of eq.(p6D. A single line suffices for the two j3 values, since the difference in Mps 
is negligible. One sees that results consistent with continuum expectations are obtained up 
to momenta of about 1.2 GeV. Even for the largest momenta deviations are less than ~8%. 



VII. SUMMARY 

In this article we presented further tests of Jb^ determinations using NRQCD heavy 
quarks and clover light quarks. We checked for scaling, studied the dependence on the 
momentum of the decaying Bg meson and also investigated cancellation of power law terms 
through 0{as/{aM)). At one value of the lattice spacing we were able to determine the 
Ds meson decay constant. Our best values for fs^ and /d^ are given in eq. (|T5|) and eq.(p2|) 
respectively. Our results are in good agreement with other lattice determinations of these 
quantities. In the case of fn^, for which experimental measurements exist, we obtain a value 
consistent with the current experimental world average as given in [pO| . 
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APPENDIX A: MATCHING COEFFICIENTS FOR A 



In Table X we collect matching coefficients for Aq, i.e. the po, pi, p2 and the mixing 
coefficient Cio defined in eg. (p!7|) and eq.(|l^). This augments numbers given in reference pO 



where, for the case of the full NRQCD action used in the present simulations, a smaller set 
of heavy quark masses was covered. The coefficient p2 in Table X includes contributions 

j(disc) 



from Jg , as explained after eq.(ll). We note that in the first article of reference |20 



different convention was used and p2 there did not include 2aMoCyio- The new convention 
adopted here is more in line with those employed in the second article of reference pO| for 
matching of A^., Vq and V^. 
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TABLES 



volume 
# configs 

csw 
a-^rup) (GeV) 



Ks 



13 = 5.7 
12^ X 24 
278 

1.567 (tad. imp.) 
1.116(12)(1^6) 
4.20(25)(5)(104) H 
0.1400 
0.1399(1) [H 



P = 6.0 
16^ X 48 
102 X 2ij.g^ 
1.479 (tad. imp. 
1.92(7) i 
2.22(11) i 
0.13755 
0.13755(13) |] 



/3 = 6.2 
24^ X 48 
144 

1.61 (nonpert.) 

2.59(l?o)(l^) PHI 
1.64(5)(l«)(t?) 11 
0.1346 
0.13466(7) m 



TABLE I. Simulation details for sections I through V. The errors on a" include statistical 
errors and those due to the chiral extrapolation of rup. aM^ is the bare b-quark mass in lattice 
units determined from the B (or Bg) meson. Kg gives the light quark mass for which results are 
presented here and Kg is the actual strange k for these f3 values based on the K meson. At /3 = 6.0 
time reversed configurations were also used and results were interpolated to Kg- 















(3 = 5.7 
















aMo 


20.0 


12.5 10.0 


8.0 6.0 5.0 


4.0 


3.5 3.15 


2.75 2.45 


2.2 


2.0 


1.7 


1.5 


1.3 


1.125 1.0 0.8 


0.6 


n 


1 


1 1 


1 1 2 


2 


2 2 


3 3 


3 


3 


4 


4 


5 


6 6 8 


10 














f3 = 6.0 
















aMo 




10.0 


7.0 




4.0 




2.7 






2.0 




1.6 




n 




1 


1 




1 




2 






2 




2 
















(3 = 6.2 
















aMo 




4.5 






2.5 










1.44 








n 




1 






3 










3 









TABLE IL Bare heavy quark masses and n- values 
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Source 


p — 0. / 


p — D.U 


P — D.Z 


statistical 


3 


3 


2 


disc. U[\aKQCD) ) 




4 


o 
z 


disc. Oiag (aAocD)) 
pert. 0(a2), 0(a2/(aAf)) 


5 


2 




11 


7 


6 


NRQCD 0((Aqcd/M)2), 0(a,AQci3/M) 


1 


1 


1 


Kg 


+4 


+4 


+4 






4 


3 


Total 


19 


11 


8 



TABLE III. Estimates of the statistical and main systematic errors, in percent, in our values 
for /b. 
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aiviQ 


„3/2 ffO),/!? 


„3/2 f (l),/lu" 


„3/2 j-(3),/lu" 


„3/2 f(4),/T7 


„3/2 f(5), /TT 






Q 

P = 


0. r 






U.oUU 


U.455(7j 


-U.lo8(4) 


-U. 229(4) 


U.1148(2U) 


-U. 1463(24) 


U.oUU 


U. 445(7) 


-U. 1448(27) 


-U. 1145(21) 


U.U583(1U) 


-U.U763(12) 


l.UUU 


U.445(oj 


-U.12U4(24) 


-U.Ud78(13) 


U.U348(d) 


-U.U47U(8) 


1.125 


U.44y(oj 


-U.11U5(23) 


-U.U535(11) 


U.U2d7(5) 


-U.U369(7) 


l.oUU 


U. 455(9] 


-U.U992(22) 


-U.U395(9) 


U.U193(4) 


-U.U274(d) 


1.5UU 


U. 452(9) 


-U.U89o(21) 


-U.U297(7) 


r\ ni A'X /o\ 

U.U141(3) 


-U.U2U5(4) 


1.7UU 


U.471(lUj 


-U.U82U(2U) 


-U.U238(d) 


U.U1U72(25) 


-U.U161(4) 


z.UUU 


U.4»4(llj 


-U.U731(19) 


-U.U179(5) 


U.UU754(19) 


-U.U117U(27) 


z.zUU 


U. 492(11) 


-U.Ufa85(18) 


-U.U153(4) 


U.UUd1d(17) 


-U.UU975(23) 


Z.45U 


U.5Uo(12) 


-U.Uo3d(18) 


-U.U13U(4) 


r\ r\r\Ar\r\/i a\ 

U.UU49U(14) 


-U.UU795(19) 


Z.75U 


U. 515(13) 


-U.U587(17) 


-U.U1U8(3) 


U.UU385(12) 


-U.UU639(15) 


o.loU 


U.52o(14) 


-U.U5oo(1d) 


-U.UU874(25) 


U.UU289(1U) 


-U.UU493(12) 


0.5UU 


U. 541(15) 


-U.U495(1d) 


-U.UU748(23) 


U.UU233(9) 


-U.UU4U5(1U) 


A r\c\r\ 

4.UUU 


U.55o(17) 


-U.U45U(1d) 


-U.UUd14(2U) 


U.UU177(7) 


-U.UU315(8) 


o.UUU 


U.5oo(21) 


-U.Uo81(15) 


-U.UU44U(1d) 


U.UU112(5) 


-U.UU2U7(6) 


D.UUU 


U.DU2(2o ) 


-U.Uooi(iO ) 


-U.UUo3i(i3 ) 


n nnnw/' ^\ 
U.UUU/ /(O) 


n nni /i fi /' /i \ 
-U.UU14d(4 ) 


Q nnn 
o.UUU 


U. 04(^0 ) 


n nofiQ/ 1 c^A 
-U.U2Do(iO j 


n nnoi q/'i nA 
-U.UU2io^iU j 


n nnn/1 q/'q\ 
U.UUU4o(^o j 


n nnno q\ 
-U.UUUou^o j 


1 n nnn 
iU.UUU 


U.0D(4 j 


n noi c/i /I ^ 
-U.U2io(^i4 j 


n nni 

-u.uui4y(^y j 


n nnnoTo/OP^A 
U.UUU2 ( 2(^20 j 


n nnnc^c^F^/'oc^^ 
-U.UUU000(^20 j 


12.5UU 


U.68(6) 


-U.U18U(14) 


-U.UU1U3(7) 


U.UUU171(19) 


-U.UUU362(2U) 


2U.UUU 


U.73(9) 


-U.U119(13) 


-U.UUU45(4) 


U.UUUU66(11) 


-U.UUU144(11) 






/3 = 


6.2 






1.44U 


U.143(2) 


-U.U195(3) 








2.5UU 


U.147(2) 


-U.U125(2) 








4.5UU 


U.152(l) 


-U.UU76(1) 









TABLE IV. Matrix elements of individual current contributions in lattice units, f^^^ is defined 
in eq. (p!3|). The (3 = 6.2 matrix elements have different normalization from those at other /?'s (see 
text). 
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aMo 


f t'pp_1pvp1 

uJ- V CI 


n* — l/a 




q* = 2/a 


q* — tt/ci 






o 

P = 


5.7 






(J.oUU 


0.ol5(Dj 


0.263(5) 




A OO A / CT \ 

0.280(5) 


0.286(5) 


U.8UL) 


0.o5o(Dj 


0.280(5) 




A O A /I / C \ 

0.304(5) 


A O 1 O / C \ 

0.312(5) 


l.UUU 


O.o8z(7j 


0.302(5) 




0.328(6) 


0.337(6) 


1.1/5 


0.399(7) 


0.314(6) 




0.341(6) 


0.351(6) 


l.oUU 


0.419(5j 






0.358(7) 


0.369(7) 


l.oUU 


0.440(9j 


U.o4o( / ) 




0.3(5(8) 


0.386(8) 


1.7UU 


n /I en/"! n^ 
0.459(10; 


0.364(8) 




n on /I /n^ 
0.394(9) 


0.406(9) 


z.UUU 


0.484(11) 


0.378(9) 




A /I 1 O / 1 A\ 

0.412(10) 


A Af\ A / t A\ 

0.424(10) 


2.2(JU 


0.500(11) 


A o o r* / a\ 

0.385(9) 




0.422(10) 


A /I O CT / 1 A\ 

0.435(10) 


Z.45U 


0.51»(lz) 


A O A/^ / 1 A\ 

0.396(10) 




A /I O CT / 1 A^ 

0.435(10) 


A /1/1A/11A 

U.44y(ilj 


Z.YOU 


r\ P'00/1 yl\ 

0.5o8(14) 


A /1A0/'1A\ 

0.408(10) 




0.450(11) 


A A r*tz ^ 1 o\ 

0.465(12) 


0.15L) 


0.560(15) 


0.421(12) 




A A f ( ' / 1 0\ 

0.466(13) 


A A OO / 1 0\ 

0.482(13) 


o.oUU 


n C7n/i '7^ 
0.579(17) 


0.432(13) 




n A Tn/i /I \ 
0.479(14) 


n A n£;/i /i \ 
0.496(14) 


4.U(JU 


0.603(19) 


A A A £i f '{ A\ 

0.446(14) 




0.496(16) 


0.514(16) 


5. DDL) 


0.64z(z4) 


0.471(18) 




A C 0£? / 1 (\\ 

0.526(19) 


A C /I / OAA 

0.546(20) 


D.UUU 


0.6/ l(z8) 


0.488(21 ) 




0.54/(23) 


0.568(24) 


Q nnn 
o.UUU 


0. /z(4) 


0.519(28 ) 




0.58(3) 


0.61(3 ) 


iU.UUU 


u. (0(0 ) 


U.00(4 ) 




U.Di(4) 


n A/i //I A 
U.D4(4 ) 


1 o p^nn 


U. / y ( ( ) 


U.O / (0 ) 




U.d4(D ) 


U.D ( (0) 


on nnn 

zu.uuu 


n Q/i /I n^ 
U.o4(iU ) 


U.Do(y ) 




n vn/ 1 n^ 
U. ( U(iU ) 


n vo/' 1 n^ 
U. ( Z(iU ) 






o 

p = 


n 
6.U 






1 £;nn 

l.bUU 


0.442(13) 


0.398(11) 




0.410(11) 


0.414(11) 


o nnn 
z.UUU 


U. 466(13) 


U.4U6(13) 




U. 421(13 ) 


U. 42/(13) 


7nn 




n A'^^ (^P.\ 

U.4Zi^iD ) 






n /I /IS/1 fiA 
U.44o^^iu j 


A nnn 
4.UUU 


n F^oi f 1 n^ 
U. 021(19 ) 


U. 442(16 ) 




U. 460(16 ) 


n A '7Q /I 

U. 4/3(16) 


7.000 


0.588(16) 


0.479(13) 




0.508(13) 


0.518(13) 


10.000 


0.615(16) 


0.498(13) 




0.528(13) 


0.540(13) 








6.2 






1.440 


0.515(8) 


0.404(4) 




0.429(8) 


0.442(8) 


2.500 


0.561(8) 


0.465(6) 




0.485(6) 


0.495(6) 


4.500 


0.602(4) 


0.505(5) 




0.527(5) 


0.536(5) 



TABLE V. Decay matrix elements f\fM in GeV^I'^ . The first column gives tree-level results, 
the remaining columns include renormalization constants using ap at three values of g*. Only 
statistical errors are shown. 
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/b, in MeV 


5.7 


217(7)(30)(24)(2)(9)(t.^^) 




217(lt^) 


6.0 


184(6)(8)(13)(2)(7)(7) 




184(19) 


6.2 


187(4)(4)(11)(2)(7)(6) 




187(16) 



TABLE VI. Quenched Jb^ with errors listed separately and added in quadrature to the right. 
The errors denote from left to right statistical (plus interpolations in Mq), discretization, pertur- 
bative, relativistic, fixing Kg Sind a~^. 





13 = 6.0 


/3 = 6.2 


volume 


163 X 48 


243 X 48 


# configs 


268 


144 


csw 


1.77 (nonpert.) 


1.61 (nonpert.) 


aMo 


2.0 


1.44 


Kg 


0.13344 


0.1346 



TABLE VIL Simulation details for section VL 





^ = 6.0 


/3 = 6.2 


aMpgrt 


2.53(8) 


1.88(5) 








(0,0,1) 


2.51(23) 


1.93(33) 


(0,1,1) 


2.83(36) 


1.77(19) 


(1,1,1) 




1.73(17) 


(0,0,2) 


2.49(11) 


1.63(16) 



TABLE VIIL aMkin extracted from pseudoscalar mesons with different momenta and com- 
parison with aMpert- Momenta are given in units of 2ti /aL. Values for the lowest momenta in 
lattice units are: ap = 0.39 {(3 = 6.0) and 0.26 (/3 = 6.2). 
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(3 = 6.0 




(3 = 6.2 




momentum 


R{p) 


i?W(p) 


R{p) 


i?W(p) 


(0,0,1) 


1.02(2) 


1.01(2) 


1.01(2) 


1.00(2) 


(0,1,1) 


1.00(5) 


0.98(5) 


1.03(3) 


1.02(3) 


(1,1,1) 






1.05(4) 


1.03(4) 


(0,0,2) 


1.08(3) 


1.05(3) 


1.10(6) 


1.07(6) 



TABLE IX. The ratios R{p) and R^^\p) of eq.(^7|) and eq.(|2^) for several momenta and /3 
values. Momenta are given in units of 27r/aL. 
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TABLE X. One- loop matching coefficients po, p\ and pi of eq.(^) and mixing coefficient 
Cio of eq.(^) for several heavy quark mass values. p2 includes 2aMo Cao in order to incorporate 
contributions from j('^'^c). The NRQCD action is that of eq.(|) - (||) with q = 1. Both the NRQCD 
and clover light quark actions are taken to be tadpole-improved with the plaquette definition of 



22 



FIGURES 



0.4 



0.3 



> 

P. 0.2 



PQ 



0.1 



0.0 



1 \ \ r 



0.0 



1 \ \ r 



<> 



J \ \ L 



<> 



o JLQCD 

o this work 



J \ \ L 



0.5 

a (GeV"^) 



1.0 



FIG. 1. versus the lattice spacing and comparison with reference g. Systematic and 

statistical errors have been added in quadrature. 
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FIG. 2. Matrix elements f'^^^VM, |/(^)| VM and l/^iil/M versus 1/Mps at /3 = 6.0. q* = 2/a 
was used to obtain |/i^f,|\/M. 
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FIG. 3. fpsVMps versus 1/Mps at /3 = 6.0 for q* = tt/o and q* = 1/a. The lines are linear 
fits to the data. Only statistical errors are shown. 
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FIG. 4. Results for fpsV^ps versus 1/Mps from different /3's. q* = 2/a is used for the 
one-loop renormalized data points. The central full line is a linear fit to the q'* = 2/a, /? = 6.0 
data. The bottom(top) dotted line comes from lower(upper) bounds on fits to q* = l/a{q* = tt/o), 
/3 = 6.0 data (see text). 
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FIG. 5. Matrix elements /(o^^M, |/(^)|VM, |/£|,|VM and + + /(5)|)/m versus 

\/Mps at /3 = 5.7. = 2/ a was used to obtain |/i^|,|\/M. 
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FIG. 6. aMkin extracted from mesons with different spatial momenta. The fancy squares to 
the left show aMpert- The horizontal lines bracket perturbative errors on aMp^rt- 
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FIG. 7. The ratio R{p) of eq.(|27|) versus momentum. The full line shows E{p)/ y/Mps- 
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